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We argue that through the Wilson lines, gauge invariance has as an effect that the hard 
functions in weighted spin-asymmetries in hadronic scattering processes are given by 
gluonic pole cross sections, rather than the usual partonic cross sections. 

1 Introduction 

The interest in the transverse polarization properties of hadrons has been increasing consid- 
erably over the last decades after the observation of large single transverse-spin asymmetries 
Am in processes such as pp— >A°^X, p^p^nX and p^p^irX, see e.g. [2]. Typically, single- 
spin asymmetries (SSA) are defined as the difference of the scattering cross sections with op- 
posite spin orientations divided by the spin averaged cross section. As such, the spin asym- 
metries are pre-eminent observablcs for the spin dependent parton distribution functions, 
such as the quark transversity function h±(x), which cannot be probed at leading order in 
inclusive deep-inelastic scattering. 

A mechanism to generate single-spin asymmetries through soft gluon interactions be- 
tween the target remnants and the initial and final state partons was first proposed in the 
context of collinear factorization [3j S] ■ Heuce, this collinear factorization formalism involves, 
apart from the usual twist- two collinear distribution functions, also twist-three collinear 
quark-gluon matrix elements. Since they contain a gauge field operator corresponding to 
a soft gluon, they are referred to as gluonic pole matrix elements. An important example 
is the Qiu-Sterman matrix element Tp(x,x) [3]. Several other mechanisms for generating 
SSA's through the effects of the intrinsic transverse momenta of the partons have also been 
proposed. For instance, in the Sivers effect the asymmetry arises in the initial state due 
to a correlation between the intrinsic transverse motion of an unpolarized quark and the 
transverse spin of its parent hadron [5]. The effect is described by a transverse momentum 
dependent (TMD) function f^ T {x : p^). However, based on time-reversal arguments it was 
expected to vanish [6]. It will, therefore, be understood that it came as a surprise when B rod- 
sky, Hwang and Schmidt in 2002 showed that leading twist spin asymmetries can anyway be 
generated in a spectator diquark model where the scattered quark undergoes an additional 
gluonic interaction with the target remnants [7] . It was soon realized that the time-reversal 
arguments that were used to predict a vanishing Sivers effect were flawed as they neglected 
the presence of the Wilson lines. The Wilson lines or gauge links IA — V exp[— igfdz-A(z)} 
are path-ordered exponentials that are required to obtain gauge-invariant definitions of par- 
ton distribution and fragmentation functions. Taking the Wilson lines properly into account 
it was found instead that time-reversal leads to the important conclusion that the Sivers 
functions in SIDIS and Drell-Yan scattering have opposite signs 
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and that finite Sivers contributions are possible after all [S]. Single-spin effects due to the 
Sivers mechanism in SIDIS have now been observed by several collaborations [9] . 

The discussion above exemplifies the important role that the Wilson lines have come 
to play. Today it is recognized that they are not just mere operators to obtain gauge in- 
variant definitions of parton distribution functions but also that they, through the Sivers 
(and Boer-Mulders [10]) effect, can themselves be regarded as sources for single-spin asym- 
metries. The Wilson lines arc also crucial ingredients in the derivation of the relation 
2M f^, X \x) = —gTF{x,x) found in ref. [llj . demonstrating that the first transverse moment 
of the Sivers function is a gluonic pole matrix clement (the same is true for the Bocr-Muldcrs 
function). The remarkable 'universality' property of the TMD Sivers functions in Eq. (TIJ is 
a prediction intrinsic to QCD which follows from the presence of the Wilson lines or, stated 
differently, from the gluonic initial and final-state interactions. Experimental verification of 
this prediction would be profound support for our understanding of the physics underlying 
the generation of spin asymmetries. It could be tested by comparing the signs of the single 
transverse-spin asymmetries in SIDIS and Drcll-Yan scattering. As stated in the previous 
paragraph, SSA due to the Sivers mechanism in SIDIS have now been observed by several 
collaborations. However, measurements of the single-spin asymmetry in Drell-Yan scatter- 
ing are lacking behind. The reason is that the lepton-antilepton pair is a relatively rare 
final state in hadron-hadron scattering compared to purely hadronic or hadron-photon final 
states. Also in those processes the gluonic initial and final-state interactions (the Wilson 
lines) leave their fingerprints, though the effects are more intricate since the hard functions 
are more complicated. However, until precise measurements for the single-spin asymmetries 
in the 'gold-plated' Drell-Yan process become available, they may also be used to test the 
formalism describing single-spin asymmetries. The role of gauge invariance in these more 
complicated scattering processes is the topic of the next section. 

2 Dijet and Photon-Jet Production in Hadron-Hadron Scattering 

In the basic hadronic processes, SIDIS, Drell-Yan scattering and e + e~-annihilation, the hard 
functions are just simple electromagnetic vertices (at tree-level). Depending on the partic- 
ular process only initial or final-state gluonic interactions contribute and, correspondingly, 
only future and past pointing Wilson lines occur. However, when going to hadronic pro- 
cesses that involve hard functions with more colored external legs, such as in hadronic dijet 
or photon-jet production, there can be both initial and final state gluonic interactions. As 
a result, the Wilson lines will also be more complicated than just the simple future and 
past pointing Wilson lines [12] ■ In fact, the final-state interactions will give rise to future 
pointing Wilson lines at each of the outgoing partons (where the representations of the color- 
matrices will depend on the particular type of parton) as is the case in SIDIS. Similarly, 
the initial-state interactions will lead to past pointing Wilson lines at the incoming partons 
as in Drcll-Yan scattering. The Wilson lines distributed over the different external partons 
of the hard function can subsequently be joined together by making a color-flow decompo- 
sition of the hard function, such that they can be connected along the color-flow lines. In 
particular this means that the full gauge link will have a different substructure for each of 
the color-flow diagrams. More importantly, it implies that each of the Feynman diagrams 
that contribute to the hard function might have a different gauge link structure [12] . For 
the TMD distribution functions this at first sight seems to complicate things considerably. 
However, for the collincar distribution functions remarkable simplifications occur. In fact, 
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for the collinear T-even parton distribution functions all process dependence of the Wilson 
lines disappears. For the T-odd distribution functions the effect of the jungle of Wilson lines 
is to relate their first transverse moments for different processes by simple color factors: 

fHl)( T \ I _ r [D(ab^cd)] f ±(i) ( s I ( o) 

J IT \ J J D(ab—tcd) G JIT W J SIDIS ' W 

and similarly for the transverse moments of the Boer-Mulders functions. The color fac- 
tors Cg are determined by the gauge links and, hence, depend on the scattering process 
ah— >cd. Moreover, they depend on the particular Feynman diagram D that contributes 
to this process [HJ [T3] ■ For that reason they are in a natural way associated to the hard 
functions in which they can, then, be absorbed. Consequently, the resulting hard functions 
will have the generic structure da[ a ^^ c d = X)r> ^c^^Ih-Lcd' which will be referred to as glu- 
onic pole cross sections and which should be contrasted to the usual partonic cross sections 
da a b—>cd. = ^ab^cd- ^ nc summations run over all Feynman diagrams D that contribute 
to the scattering process ab^cd and the da^ are their pQCD expressions. The gluonic 
pole cross section has a bracketed subscript indicating which parton is associated with the 
gluonic pole matrix element. They appear as the hard functions whenever gluonic pole ma- 
trix elements contribute, such as the first moments of the Sivers or Boer-Mulders functions. 
This is typically the case in weighted azimuthal spin asymmetries. Hence, the effect of the 
gauge links (or, equivalently, of the initial and final state interactions) for weighted SSA's is 
that the hard functions of these observables will be given by the (manifestly gauge invariant) 
gluonic pole cross sections [13], rather than the usual partonic cross sections as might be 
expected in a 'generalized parton model approach'. Several other recent theoretical studies 
of single-spin asymmetries seem to point in a similar direction |14[ I15[ 116] . though more 
research is required to convey the exact connection between the different formalisms. 

The effects of the gluonic initial and final-state interactions for the fully TMD treatment 
of these processes is not so clear-cut. In ref. [T^] a TMD factorization formula was proposed 
for the quark-Sivers contribution to the SSA in dijet production in proton-proton scattering. 
This formula involves the gluonic pole cross sections found in refs. |13| as hard functions, 
but folded with the TMD distribution functions as measured in SIDIS (i.e. with a future 
pointing Wilson line in their definitions). On the other hand, in refs. [T2][T3] it was observed 
that complicated Wilson line structures occur in the TMD distribution (and fragmentation) 
functions in such processes. The work of [HI [13] and [16] could be related by extending 
the relation in Eq. ([2]) between the first moments of the Sivers functions to the full TMD 
functions, in analogy to the relation between the Sivers functions in SIDIS and Drell-Yan 
scattering as given in Eq. ([1]). However, the jungle of Wilson lines found in refs. [12] [13] in 
concurrence with a model calculation led the authors of ref. [T7] to conclude that a TMD 
factorization formula for spin asymmetries in processes such as dijet production in proton- 
proton scattering is not possible with universal distribution functions. It is also asserted 
that a proof of TMD factorization for such processes will be essentially different from the 
existing proofs for SIDIS and Drell-Yan scattering and that it will possibly involve 'effective' 
distribution functions [15| . Future research will have to clarify what the relation is between 
the results of refs. [HI Q2] HH E] . In particular, soft factors might play an important role 
in such a comparative study. 

Using the azimuthal imbalance of the outgoing jet pair (p^p—>jjX) or photon-jet pair 
(p* p—ryjX) weighted single-spin asymmetries have been constructed that involve the gluonic 
pole cross sections at leading twist [T3l [HI • Predictions for the quark-Sivers contribution 
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to these spin-asymmetries in dijet production have been presented [TJ [TS], making use of 
recently obtained parametrizations |20j for the quark-Sivers first-moments (Qiu-Sterman 
matrix elements) obtained by fitting to An data taken at E704, BRAHMS and STAR, and 
for the quark-Sivers half-moments obtained by fitting to SSA's measured at HERMES. First 
experimental measurements for this process have been performed at RHIC in 2006 [2Tj . The 
results seem to indicate that the full QCD treatment, which makes use of gluonic pole cross 
sections to account for the competing effects of the gluonic initial and final-state interactions, 
is consistent with the data, while the 'generalized parton model approach', which only uses 
partonic cross sections, overestimates the measurements |21| . However, the experimental 
results are at the same time consistent with vanishing SSA's in dijet production and more 
measurements to increase statistics and disentangle flavor dependence in concurrence with 
more elaborate theoretical studies would therefore be welcome. 

Predictions for the single-spin asymmetry in the similar photon-jet production process 
have been made in kinematical regions accessible at RHIC. For these predictions use was 
made of parametrizations [22] for the quark-Sivers function obtained by fitting to HERMES 
and COMPASS data. The main result is that the effect of the gluonic initial and final-state 
interactions is to flip the sign of the spin-asymmetry in the full QCD treatment compared to 
the generalized parton model result jT] 119) . This can be understood intuitively by observing 
that the dominant partonic channel is qg—>"/q scattering, where in the large-iV limit the 
color flows from the initial quark and through the gluon back into the initial state. There 
the initial-state interactions lead to the sign flip of the spin asymmetry, as is also the case in 
Drell-Yan scattering. However, in the qg—>"fq process the initial-state interactions are gluon 
self-couplings and, hence, the experimental measurement of the predicted sign flip in this 
case would test the very non-abelian nature of QCD. It is also shown that the Boer-Mulders 
and gluon-Sivers effects are too small to affect this conclusion. Therefore, by measuring the 
sign of the asymmetry one has a test of the QCD formalism to describe the generation of 
single-spin asymmetries that is equally legitimate as the measurement of the SSA in Drell- 
Yan scattering. 
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